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We present an analytic study for subdiffusive escape of overdamped particles out of a cusp-shaped 
parabolic potential well which are driven by thermal, fractional Gaussian noise with a l/u> 1 ~ a power 
spectrum. This long-standing challenge becomes mathematically tractable by use of a generalized 
Langevin dynamics via its corresponding non-Markovian, time-convolutionless master equation: We 
find that the escape is governed asymptotically by a power law whose exponent depends exponentially 
on the ratio of barrier height and temperature. This result is in distinct contrast to a description with 
a corresponding subdiffusive fractional Fokker-Planck approach; thus providing experimentalists an 
amenable testbed to differentiate between the two escape scenarios. 

PACS numbers: 05.40.-a, 82.20.Uv, 87.16.Uv 



The theme of anomalous sub-diffusion and rate kinet- 
ics continuous to flourish over the last years. This topic is 
driven by the availability of a wealth of intriguing exper- 
imental data, ranging from anomalous diffusion in amor- 
phous materials, quantum dots, protein dynamics, actin 
networks, and biological cells [l], [1, H, 13, H B 0> B H> 03 
[IlT ] . Suitable theoretical descriptions derive from contin- 
uous time random walks (CTRW) 0, 0, the CTRW- 



based fractional Fokker-Planck (FFP)-approach |13|, 14 1, 
or the generalized Langevin equation (GLE) [15|, KL6 |. 
The GLE-subdiffusion implies power-law-correlated ther- 
mal forces (or fractional Gaussian noise (fGn) 17]) pos- 
sessing infinite memory with a \/uj l ~ a power spectrum 
@, lle| Such random forces emerge when coupling the 
system to sub-Ohmic thermal baths with spectral densi- 
ties J(cj) oc co a , < a < 1 [HI]. 

Recent experiments on anomalous conformational sub- 
diffusion in electron-transferring proteins have been suc- 
cessfully modeled within a FFP equation [f| . Soon after, 
however, the fGn-Langevin approach has been shown to 
describe the experimental data even more convincingly 
Q. Both approaches are consistent with molecular dy- 
namics simulations (l9j . In addition, both schemes are 
consistent with the laws of thermodynamics. Fundamen- 
tal differences manifest themselves, however, when one 
considers the escape dynamics. The description of fGn- 
driven escape presents a long-standing, timely problem. 
This is so because this complexity, in contrast to the 
Kramers escape dynamics [15j . then generally no longer 
allows for a detailed, even approximate solution. Re- 
cent attempts to solve this challenge can be found in 
Refs. [2(3], although not yielding a proper solution. Thus, 
the issue of 1 /f-noise driven escape remains an intriguing 
open problem which continues to haunt the literature. 

Our analytic work deals with the unique, but exactly 
solvable challenge of an escape driven by 1 /f-noise out of 
a cusp-shaped parabolic potential. In doing so we demon- 
strate that the resulting escape dynamics is scale-free, i.e. 
it is governed by a power law. This in turn invalidates a 
rate description. Moreover, the escape dynamics within 
the fGn-Langevin description is exponentially sensitive to 
temperature. This result is in marked contrast to the de- 



scription within a FFP-dynamics which instead yields a 
temperature-independent power law exponent. Thus, in 
contrast to the finite mean first passage time (MFPT) 
result in [2(| our main result exhibits an infinite MFPT, 
being consistent with a strict subdiffusive escape dynam- 
ics. A rate description emerges only when invoking a 
physically plausible low frequency regularization of the 
noise spectrum. 

GLE-approach. We start out from the GLE for a par- 
ticle of mass m moving in the potential V(x) (l5| : i.e., 



mx + / rj(t - t')x(t')dt' + V'{x) = £(t). 



(1) 



The autocorrelation function (£(i)£(i')) of thermal Gaus- 
sian noise £(t) and the frictional kernel rj{t) are related 
by the usual fluctuation-dissipation relation [Til ]: 



(Z(tW))=kBTr ) (\t-t'\). 



(2) 



In the following we consider the overdamped limit with 
m — * 0; i.e. the velocity is thermally relaxed at each 
instant of time. Moreover, we assume that the particles 
are initially localized in a metastable parabolic well at xo, 
cf. the inset in Fig. ([T]). The starting probability den- 
sity then is P(x, t = 0) = 5(x — xq). The corresponding 
non-Markovian master equation for P(x,t) for this GLE 
is generally not known: For arbitrary physical memory- 
friction rj{t) this task is known for two cases only; namely 
(i) a linear potential V(x) = —Fq x, including the case 
of free diffusion, i.e. Fq = and (ii) a parabolic poten- 
tial V(x) = kx 2 /2. The procedure to obtain the master 
equation is well known: It is so lely rooted in the Gaus- 
sian nature of x(t) [II, EI El El El ■ The result is a 
time-convolutionless master equation for P(x, t) obeying 
a Fokker-Planck form with a time-dependent diffusion 
coefficient D(t) [HHH, reading 



dP{x,t) _ 8_ 
dt ~ [) dx 



ox 



(3) 



Here, (3 = l/(ksT). Notably, D(t) does not depend on 
xq, but is dependent on V{x) and memory friction rj(t). 
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For a quadratic potential, it can be expressed via the 
relaxation function 9(t) of position fluctuations as [22LI24T] 



D(t) 



dt 



ln0(*) , 



(4) 



where It — y/ksT/n is the length scale of thermal fluc- 
tuations. The Laplace-transform of 0(f) is related to the 
memory friction 77(f) by 6(s) = fj(s)/[K + sfj(s)] . 

We next use a power-law friction kernel 77(f), reading 



v(t) = 



r(i 



< a < 1 



(5) 



This friction 77(f) yields an anomalous, free (V(x) = 0) 
subdiffusion with (5x 2 (t)) = 2K a t a /T(l + a), where the 
anomalous diffusion coefficient K a = ksT/rja obeys a 
generalized Einstein relation. For a parabolic potential 
this yields the relaxation function 



6(t) = E a [-(t/r D ) a ] 



(6) 



where td — (^a/^) 1 j and E a (z) is the Mittag-Lcfllcr 
function, i.e., E a (z) — J2^=o z n /T(an+ 1) [26[. It corre- 
sponds to the Cole-Cole model of glassy dielectric media 
[271 ] . whereas the limit a — ► 1 corresponds to an expo- 
nential relaxation with E±(x) = exp(x). 

The thermal fGn £(*) is the time derivative of frac- 
tional Brownian motion (fBm) 17] with a power spec- 
trum, S((w) = 2fc B Tr/ Q sin(7ra/2)/w 1_Q . For this ther- 
mal fGn the CLE in |T]) with m — can formally 
identically be recast as the "fractional" Langevin equa- 
tion, i.e., r) a D%x(t) + V'(x) = £(*), wherein D°x(t) = 
(l/r(l - a)) J 'd*'(* - t')- a Mf) is the operator of the 
fractional Caputo derivative [261 ] . 

FFP- approach. Alternatively, if instead of the fGn 
dwelling in a potential in |T]), ([2]), @we use a modeling in 
terms of an overdamped, fractional Fokker-Planck equa- 
tion description [ID, [3] , the probability density obeys 



D°P(x,t) = K c 



d_ 

dx 



-pv {x) d 0v[x) \ _ 

ox J 



This result derives from an underlying continuous time 
random walk description of subdiffusion [lj. It as well 
has an associated Langevin equation in a random oper- 
ational time *(t) [28| which, however, is profoundly dif- 
ferent from the GLE. 

Subdiffusive dynamics dwelling in a parabolic poten- 
tial. The time-convolutionless master equation ([3]) of 
the GLE in Q, ©, dSJ) can be solved exactly for 
a parabolic potential [29(: We first transform P(x,t) 
as P(x,t) — exp(— (3V (x) / 2)W (x , t) and separate the 
variables, W(x, *) = Y"(x) ( I > (t). For the coordinate- 
dependent part this yields a spectral representation, 
reading 

Y'n{x) + ~ k (l - ^ra 2 + 2X n /((3 K )\ Y n (x) = 0, (8) 



where A„ and Y n (x) are the corresponding spectral eigen- 
values and eigenfunctions. The functions <&„(*) obey: 

$n(*) = -A„£>(*)$„(*). (9) 

By use of (U]) the exact solutions of ^ read 

<!>„(*) = [0(*)] s », (10) 

where s n := l T X n . These findings yield for the explicit 
solution for the probability density P(x,t) the result, 



P(x,t) = exp(-f3Kx 2 /4)Y,CnY n (x)[e(t)Y 



(11) 



where the expansion coefficients c„ are determined from 
the initial probability density P(x,t = 0). The spec- 
trum reads s„ = n, n — 0, 1,2, .... Moreover, the func- 
tions Y n (x) are given in terms of Hermite functions [l5| . 
The dynamics of the probability evolution is thus ruled 
by the relaxation function 0(f) in Eq. ((6]). Remark- 
ably, the relaxation of the mean value (x(t)) follows 
precisely to the same law as in the case of a FFP de- 
scription p^| . This is surprising because the general 
solution P(x, t) differs markedly from that of the FFP 
equation in [lj]. The solution in the latter case is ob- 
tained from (fTTj) by substituting therein E a [— n(f/T£>) a ] 
for (E a [-(t/T D ) a ]) n within 9 n (t). This follows from the 
fact that for the FFP in Q, equation Q is replaced by 
D?$ n (t) = -\ n K a $ n (t) M- 

Because the integral over 0(f) in Eq. ^ is not finite, 
the statistical mean of random escape times is subject 
to a divergence, thus invalidating a rate description. Be- 
cause inter-well transitions can always be broken up into 
the two steps, i.e. (i) reaching the barrier region from 
well bottom and (ii) a barrier (re)-crossing to an adja- 
cent well, this implies as well a diverging MFPT when 
considering more generic situations as the stylized one 
addressed next. 

Escape out of parabolic cusp potential. Let us impose 
next an infinitely sharp potential cutoff at x a — L 3> It, 
see the inset in Fig. [TJ This sharp cut-off is identical to an 
absorbing boundary condition, satisfying P(x, t) = for 
x > L. The Gaussian approximation for the GLE there- 
fore remains valid inside the parabolic cusp potential [30l ] . 
The solution of eq. ([3]) for the corresponding boundary 
value problem now reads anew (for — oo < x < L): 



Y n (x) = U -s n 



-x/l T 



(12) 



where U(u,x) denotes the parabolic cylinder function. 
The spectrum for this case reads A„ = s n /l T , being de- 
termined by the solutions of transcendental equation 



r (s n -^-L/i T ) =o. 



(13) 



For L — ► oo, s n approaches again n. For large L 3> It the 
decay of the survival probability inside the well Pgp (*) = 
(x,t)dx is ruled by the lowest eigenvalue; i.e., 



/-oo ^cuspi 



P S p(*)«[0(*)]*°. 



(14) 



This constitutes our first central result. The value so is 
given by the numerical solution of Eq. (flU)) . Remark- 
ably, it is well approximated by the inverse of the prop- 
erly scaled mean first passage time of the correspon ding , 
memoryless Markovian problem yielding, cf. in Ref. [la ]: 



h- 1 - I- 2 

a — Irp 



dy 



dxexp(/3[V(y)-V(x)}). (15) 



This yields Sq = F 



where Vo 



kL 2 /2 is the 



barrier height and l/F(z) = ^/tt e y [1 + erf(y)]dy. 
For example, for L — 21t the exact value of sq is 

s (oxact) = q 0g727 whik yieMs s (approx.) = q ggggg 

The difference is already less than 1.5% and rapidly di- 
minishes with increasing L. As a main trend, sq de- 
creases approximatively exponentially oc exp(— Vo/fc^T), 
thus displaying a typical Arrhenius dependence. 

Within the approximation of (QT]) the MFPT of the 
non-Markovian escape dynamics is given by: 
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FIG. 1: Comparison between the escape dynamics de- 
scribed by the fractional Fokker-Planck equation and the 
overdamped, non-Markovian GLE-dynamics in HJ-([5]| in a 
parabolic cusp potential, depicted with the inset. 



[T) I P SP (t)dt 



[6(t)] s °dt , 



(16) 



being indeed very distinct from the Markovian case. For 
the fGn-GLE model, denoted in the following by Pgp LE , 
we find from JTU) P^ P LE f) w (E a [-(t/T D ) a ]) s ° , and thus 
the MFPT diverges, i.e. (r) = oo. 

Likewise, the high-barrier solution of the FFP equation 
in © is given by P S F P FP (0 » E a [- S() {t/T D )% yielding 
again no finite value for the MFPT. The asymptotic long- 
time behaviors differ distinctly in these two models: 



pGLE 



(t) 



1 (1E.Y 
s r(l-a) V t J 



(17) 
(18) 



where td = iVa/n) 1 ^ "- In particular, for the fGn-GLE 
model, the power law exponent squ depends exponen- 
tially on the barrier height and the (inverse) temperature. 
In contrast, for the FFP-theory this power law exponent 

just equals the subdiffusive power law exponent a. 

Fig. [1] depicts a comparison between two es- 
cape dynamics for a = 1/2 Q, where 9 if) = 
exp(i/TD)erfc( v /</r D ), and for (3V = 2. The FFP- 
escape dynamics overall proceeds faster. The initial 
kinetic stages are identical. A detection of a power- 
law escape that is exponentially sensitive to tempera- 
ture would corroborate the GLE based approach; while 
a temperature- independent power law decay would favor 
the FFP-approach. 

Role of memory cut-off. Both considered theoreti- 
cal models contain a physical drawback: Random forces 
obeying a true l/w 1- " feature in the power spectrum 
are not physical; i.e. a low-frequency regularization 
must always emerge on physical grounds, (implying that 
S(w = 0) is finite) [25]. To account for this physical re- 
quirement we introduce an exponential cutoff exp(— uj c f) 




w 2 10° io 2 io 4 io 6 io 8 io 10 io 12 io 14 io 16 

time, in units of x D 

FIG. 2: Survival probability for a modified fGn-GLE model 
using a regularization. An exponential relaxation tail emerges 
for t > LUc 1 . The chosen parameters are: td ~ 1 ps, lo c ~ 
1 hour^ 1 and a — 1/2. 



with a small frequency lu c for the memory kernel ([5]), 
yielding S(lj = 0) = 2kBTr/ a /ujl~ a . The correspond- 
ing relaxation function 9 if) now exhibits an exponential 
decay 9(f) oc exp(— ui c t), for t > l/u) c . The memory ker- 
nel rj(t) becomes integrable so that the MFPT (r) exists. 
As a consequence, a non-Markovian rate description now 
becomes valid (l5j . In practice, however, a feasible rate 
description fails whenever the main part of the escape dy- 
namics occurs within a distinct non-exponential, power 
law regime which extends over many temporal decades. 
This feature is elucidated with Fig. [5J A valid rate de- 
scription, although with an extremely small rate is re- 
stored by either lowering the temperature, or likewise, 
by increasing the barrier height. The top curve (high- 
est barrier case) in Fig. [21 depicts this trend. Our re- 
sults corroborate also with the numerical simulations of 
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bistable dynamics [3l| . where a numerical cut-off is in- 
trinsically present. Using a memory cut-off within the 
CTRW description for the FFP in Eq. (|7|) does result as 
well in an asymptotically finite rate. The intermediate 
power law will exhibit, however, also no distinct temper- 
ature dependence, being again in a clear contrast with 
the subdiffusive GLE description. 

In conclusion, we put forward an analytical treatment 
of the survival probability for the non-Markovian escape 
from a cusp-shaped well when anomalous subdiffusion 
is acting. Then, the MFPT diverges which in turn in- 



validates a rate description. The sensible physical re- 
quirement of a low-frequency regularization enables one 
to restore a rate theory description that is valid for suf- 
ficiently high barriers, or very low temperatures. The 
single-molecular enzyme kinetics [6|, [7| might present a 
suitable candidate to validate experimentally the intrigu- 
ing crossover between an exponential and a power law 
kinetic regime which crucially depends on temperature. 
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